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In epilepsy it has been challenging to detect early changes in brain activity that occurs prior to seizure onset 
and to map their origin and evolution for possible intervention. Here we demonstrate using a rat model of 
generalized epilepsy that diffuse optical tomography (DOT) provides a unique functional neuroimaging 
modality for noninvasively and continuously tracking such brain activities with high spatiotemporal 
resolution. We detected early hemodynamic responses with heterogeneous patterns, along with intracranial 
electroencephalogram gamma power changes, several minutes preceding the electroencephalographic 
seizure onset, supporting the presence of a "pre-seizure" state. We also observed the decoupling between 
local hemodynamic and neural activities. We found widespread hemodynamic changes evolving from local 
regions of the bilateral cortex and thalamus to the entire brain, indicating that the onset of generahzed 
seizures may originate locally rather than diffusely. Together, these findings suggest DOT represents a 
powerful tool for mapping early seizure onset and propagation pathways. 

Epilepsy is one of the most common and devastating neurologic diseases affecting over 2.5 M Americans. 
Seizure control with antiepUeptic drugs (AEDs) has been proven to decrease morbidity and mortality and 
remains the only treatment option for the majority of patients. Nevertheless, despite progress in antiepi- 
leptic drug development, at least 30% of all epilepsy cases remain resistant to current therapeutics''^. For patients 
without complete seizure control, the sudden, unforeseen occurrence represents one of the most debilitating 
aspects of epUepsy. Amongst the types of epUeptic seizures, generalized tonic-clonic seizures are the most 
dangerous form, which has the potential to cause severe injuries or even death'''*. However, by providing 
appropriate and timely prevention, the risk factor can be significantly reduced\ 

In the absence of completely controlling a patient's epilepsy, predicting seizures is an important goal of clinical 
management and treatment. A question of fundamental importance to seizure prediction is to better understand 
the brain dynamics during the transition from interictal to the seizure state. A promising and widely studied 
candidate for seizure precursor is high frequency oscillations (HFOs). Recent studies have shown that HFOs 
occur in advance of impending seizures indicating its strong relationship with the seizure onset zone'"''. 

Most seizure prediction methods are based on mathematically extracting specific features from multiscale 
intracranial electroencephalogram (lEEG) recordings'*. Although studies have focused on prospectively testing 
these methods, no study to date has yet confirmed the ability of any method to predict seizures better than 
random'^. Invasive intracranial electrodes are required to detect oscillations generated in very localized or deep 
brain regions. Owing to the lEEG's insufficient spatial sampling, a barrier to obtaining better seizure predictive 
rates has been the inconsistency in identifying the earliest changes prior to seizure onset. A deeper understanding 
of the regional interactions at the beginning and during the evolution of a seizure may help identify the earliest 
seizure susceptible brain regions'" 

To advance our knowledge about seizure initiation, functional brain imaging techniques like single-photon 
emission computed tomography (SPECT) and functional magnetic resonance imaging (fMRI) are used to answer 
questions regarding whether so-called 'preictal states' can be distinguished from the ictal period and the func- 
tional networks in the brain responsible for seizure generation. With SPECT's ability to capture a "snapshot" of 
cerebral blood flow, studies in"''* observed cerebral blood flow (CBF) increases a few seconds to several minutes 
preceding clinical electroencephalographic (EEG) seizure onset in temporal lobe partial seizures. With fMRI's 
high spatial resolution and its combination with EEG, fMRI studies recently reported that blood-oxygen-level- 
dependent (BOLD) changes occur minutes prior to seizure onset, either in animals or in humans'*''^'"'. However, 
due to relatively low temporal resolution and lack of portability, such techniques are not suitable for rigorous 
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statistical seizure prediction analyses, which often require continu- 
ous longtime monitoring of the period prior to, during, and after 
seizures, as well as the continuous interictal state". 

Taking advantage of the near-infrared (NIR) window in tissue, 
optical imaging brings new opportunities to brain studies. Changes 
in optical spectroscopic properties of brain tissue are often correlated 
with changes in neuronal activity". Activity evoked optical spectro- 
scopic changes measured in the brain are most likely generated by 
changes in cerebral hemodynamics. Therefore, such optical prop- 
erties changes are also the hallmarks of seizure generation and the 
spread of ictal neural activity. With the ability to provide excellent 
brain surface maps, intrinsic optical signal (lOS) methodologies have 
been used in animal models to study seizure foci localization, pro- 
pagation, and neurovascular coupling"* ''^. It was also found in an 
intraoperative procedure that local change of intrinsic signal pre- 
ceded the seizure onset by approximately 20 seconds^". Functional 
near-infrared spectroscopy (fNIRS) has also been applied to invest- 
igate epileptic seizures""'^^ and other brain disorders^"*, especially in 
daily clinical use due to its non-invasive and portable merits. A major 
limitation of lOS and fNIRS is that each fails to provide much depth 
information. However, depth information can be obtained by using 
tomographic reconstruction methods such as diffuse optical tomo- 
graphy (DOT)^'*'^^. Here, multispectral NIR light is sequentially deliv- 
ered to the tissue at multiple locations and the diffusely scattered light 
are measured at multiple positions along the scalp, which are used to 
reconstruct the spatial distribution of tissue absorption and scatter- 
ing coefficients at each wavelength through a light propagation 
model in tissue. The images of absorption spectra are then used to 
derive the images of functional parameters such as oxy- and deoxy- 
hemoglobin^*". Indeed, DOT offers sub millimeter resolution within 
centimeter thick tissue^^. In the area of neuroimaging, DOT has been 
applied to study hemodynamic response to motor and visual stimu- 
lations in several cortical areas^*'^'^, and to quantify cerebral per- 
fusion™ and blood flow", both in humans and small animals. 
These studies show that DOT is capable of detecting hemodynamic 
responses spatially equivalent to fMRI. 

The primary goal of this study was to identify hemodynamics prior 
to and during seizures in a rat model of generalized epilepsy. Our fast 
three-dimensional (3-D) DOT system, combined with simultaneous 
lEEG recordings, allowed us to track both changes in hemodynamics 
and local field potentials (multiunit activity) throughout the brain 
during Pentylenetetrazol (PTZ)-induced generalized tonic-clonic 
seizures in rats. We found that hemodynamic changes precede seiz- 
ure onset by ~ 1 to 3 minutes, and that seizure zones initiated from a 
regional onset and then distributed to the entire brain. We also 
observed the decoupling between local field potentials and hemody- 
namic changes, accompanied by concurrent scattering changes. For 
the first time, the whole process of generalized tonic-clonic seizures 
was recorded in vivo, including the interictal state, the preictal trans- 
ition state, and the ictal onset state. Our findings suggest that a 
pre-seizure state can be detected with optical imaging and the "gen- 
eralized" seizures originate from a regional (focal) onset zone. 

Results 

Hemodynamic responses elicited during generalized PTZ-lnduced 
seizures. Generalized tonic-clonic seizures and corresponding hemo- 
dynamic changes were observed in 9/9 seizure rats elicited by i.p. 
PTZ. Time courses of the EEG gamma power and hemodynamic 
changes for all nine experiments are shown in supplementary Fig. Si. 
A representative (rat #1) co-located neural and hemodynamic 
response from left primary motor (Ml) is shown in Fig. 1. Time 
series data of local concentration of oxy- [Hb02] and deoxy- 
hemoglobin [HbR] were extracted from the same location as the 
lEEG recording. Continuous spike and wave discharges occurred 
abruptly within 4 minutes after i.p PTZ (Fig. la). Power spectrum 
density changes can also be seen on the time-frequency analysis 



(Fig. lb). Significant hemodynamic changes (t-test, p < 0.05) are 
displayed in Fig. Ic. During the resting state, [Hb02] and [Hb^] 
showed no significant changes, only with little fluctuations around 
their resting state values of —30 |iM and —10 |.iM, respectively. Chro- 
mophore concentrations were calculated from the reconstructed 
absorption coefficients distribution from three wavelengths. After 
PTZ injection, a 15 |J,M decrease of [Hb02], and a corresponding 
3 |xM increase of [Hb^] was observed. The hemodynamic responses 
were comprised of a large change of [Hb02] with a relatively smaller 
change of [Hb^]. These findings are consistent with other studies of 
epilepsy using optical methods like NIRS'^ diffuse optical imaging 
(DOI)'''', where [Hb02] and the concentration of total hemoglobin 
[HbT] followed similar trends accompanying by a small variation of 
[Hbij]. Interestingly, lEEG power increase (Fig. lb) in the fast gamma 
frequency range (80 — 120 Hz) was observed before seizure onset 
(Pearson Correlation: 0.89) with the hemodynamic change. Similar 
phenomenon has been reported before by comparing optical intrinsic 
signals and local field potentials (LFPs)'''. In the age-matched sham 
saline control experiments, no detectable lEEG spikes or significant 
hemodynamic changes were observed (Fig. Id-f). 

Hemodynamic changes precede lEEG spikes. We also found that 
hemodynamic changes began before seizure onset in 9/9 rats, 
although the pre-seizure time (lead time) for observed changes 
varied for each seizure. A statistics of such neural and hemodyna- 
mic changes for all nine rats experiments is listed in the supple- 
mentary Table SI. The lead time represents the interval between 
the moment when the first significant hemodynamic changes (t- 
test, p < 0.05) appear in any voxel and the electrographic onset 
time. A representative result from rat #2 is shown in Fig. 2. 
Locations of two regions of interests (ROIs), one from the motor 
cortex and the other from temple lobe are shown in a coronal slice 
(Fig. 2a left). By tracking the change of [HbT], whereas a significant 
hemodynamic decrease (t-test, p < 0.05) in ROIl (Fig. 2c, top) and 
increase in R0I2 (Fig. 2c, bottom) occurred about 20 seconds after 
PTZ injection, the lEEG (recorded at ROIl) seizure onset began 1 
minute after PTZ injection (Fig. 2b). The gamma power, having a 
tight relationship with [HbT] (Pearson correlation, 0.85), started to 
change almost at the same time as [HbT] (Fig. 2b). In Fig. 2d, the 
activation maps of [HbT], superimposed on the MRI template, 
shows the evolution of hemodynamic distribution in the whole 
brain. We use cool or warm color to specify [HbT] decreasing or 
[HbT] increasing. [HbT] started to change in the preictal period and 
remained relatively stable after tonic-clonic seizure onset. In Fig. 2e, 
we compare the results from all experiments (n = 9) of the time 
lengths between PTZ injection and occurrence of electroencepha- 
lographic seizure onset, significant gamma power changes and 
significant hemodynamic changes. On average, the hemodynamic 
change, similar to the gamma power change, leads lEEG seizure 
onset —80 seconds. 

Decoupling of lEEG and hemodynamic changes. Both neural (local 
field potentials) and hemodynamic activity (HbT, Hb02, and Hb^) 
from bilateral primary motor cortex (Ml), bilateral somatosensory 
cortex trunk region (SlTr), bilateral hind limb primary somato- 
sensory cortex (SIHL), and bilateral anterior thalamic regions, 
were studied using simultaneous lEEG and DOT (Fig. 3a). Results 
of seizure experiment from rat #2 are shown here. We note that a 
local dissociation between hemodynamic changes and local field 
potentials (left Ml, SlTr, Thalamus and right Ml) during the 
preictal period was observed, suggesting a lag between measurable 
local neural activity (gamma power) and local hemodynamics 
relative to seizure onset. Moreover, changes of reconstructed local 
scattering coefficients, as much as 0.3 mm" ' compared to the resting 
state value of —1.3 mm"', occurred concurrently with the hemo- 
dynamic changes. These findings further suggest that changes in 
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Figure 1 | Neural and hemodynamic responses elicited during generalized PTZ-induced seizures, (a), (b) and (c) are results from seizure experiments of 
rat #1 with PTZ injection sharing a same time line; (d), (e) and (f) are results from control experiments with saline injection sharing a same time line. The 
time point of PTZ or saline injection is indicated by the red dots. EEG and hemodynamic traces are picked up at the same location in left Ml. (a) Seizure 
spikes were observed from lEEG recordings after PTZ injection, (b) Obvious power spectral density changes EEG power increase in the fast gamma 
frequency range (80 ~ 120 Hz) were seen in lEEG time-frequency analysis during seizures, (c) Hemodynamic responses of [HbT], [Hb02] and [Hbn] 
from left Ml induced by PTZ. (d) lEEG recordings in the control experiment, (e) lEEG time-frequency analysis for control experiment, (f) No obvious 
hemodynamic changes were observed in the control experiments. 



hemodynamic activity alter tissue scattering properties, not readily 
visible with concurrent lEEG. 

We also observed early gamma power increase during the preictal 
period from most lEEG electrodes. This may indicate an early 
buildup of synchrony as a preparation for the generalized seizure 
onset. During seizures, we can see that seizure spikes with larger 
intensities (left side) tend to accompany with larger hemodynamic 
changes, suggesting a correspondence between local neural and 
hemodynamic activities as seizures progress into ictal stage. 
Besides, despite the existence of seizure spikes, the hemodynamic 
changes in SIHL (left and right) and SlTr (right) have neither obvi- 
ous increasing nor decreasing trend in such areas. However, large 
fluctuations in these areas are still notable, indicating that the hemo- 
dynamic activities are not as uniformly distributed as neural activ- 
ities during the "generalized" seizures. 

Evolution and distribution of seizures. Amongst all nine rats, we 
selected three experiments from rat #2, #3 and #4 with increasing 
seizure intensities to show the overall seizure evolutions through the 
whole brain. Fig. 4a gives the positions of five coronal slices selected 
for activation maps. The coordinates of these slices range from 
— 5 mm to 5 mm relative to bregma. lEEG recordings and 
statistical activation maps that depict [HbT] increases (warm 
colors) and decreases (cool colors) compared with the baseline 
DOT images (t-test, p < 0.05), are shown in Fig. 4b-d. 

Seizures activities shown in Fig. 4c were more robust than that 
depicted in Fig. 4b in term of amplitude and frequency of spikes. 
Accordingly, larger hemodynamic changes, positive or negative acti- 
vations were observed in Fig. 4c compared to changes in Fig. 4b. 
Amongst seizures Fig. 4d shows a high grade seizure in which 



continuous and high amplitude polyspikes were seen with concom- 
itant behavioral changes twitching limbs and stretching its body even 
under anesthesia. Furthermore, at the peaks during the preictal per- 
iod, the activation regions in Fig. 4b,c remained relatively stable 
during the ictal period. Conversely, in Fig. 4d, as lEEG spike pattern 
changed strongly during ictal period, the activation regions changed 
accordingly. We observed large hypoperfusion area in Fig. 4d, pos- 
sibly due to the excessive energy consummation of the electrical 
activities during the fierce seizures onset. 

Discussion 

Multispectral DOT was used to investigate seizure dynamics in a rat 
PTZ seizure model. Changes in local and global hemodynamics and 
local field potentials were identified during seizure initiation and 
propagation. Early hemodynamic changes with concurrent gamma 
power increase were found several minutes prior to seizure onset 
supporting the presence of a preictal state. The early changes in 
optical scattering property offer a possible explanation for the 
decoupling of local field potentials from early hemodynamic changes 
during the preictal stage. We also found that generalized seizures 
initiate with focal hemodynamic changes and subsequently evolve 
to widespread changes throughout the brain. 

There is general agreement that despite pharmacological, neuro- 
modulation, dietary, and surgical advances in the treatment of 
chronic epilepsy, seizures cannot be controlled in as many as 25% 
patients, and therefore there is a need for new therapeutic 
approaches'". At the same time, there is growing awareness that the 
development of new therapies has slowed, and to move toward new 
and more effective therapies, novel approaches to therapy discovery 
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Figure 2 | Hemodynamic changes precede EEG seizure onset, (a) Locations of two ROIs, one from the motor cortex and the other from temple lobe, 
shown on a coronal slice (right). The sagittal plane (left) gives the position of the coronal slice, (b) EEG recording from electrode placed at ROIl and its 
gamma power change. EEG seizure onset began about 3 minute after PTZ injection. EEG gamma power changes started much earlier than EEG spikes, (c) 
[HbT] changes (blue) with fitted curves (red) from two typical ROIs. Significant hemodynamic changes started 20 s after PTZ injection, (d) Activation 
map of HbT showing the evolution of hemodynamic distribution, superimposed on the MRI template. Negative activation (cool colors) means 
decreasing [HbT] , while positive activation (warm colors) stands for increasing [HbT] . Blue lines specify the time period corresponding to each activation 
map. (b) and (c) share a same time line, (b) ~ (d) are results from rat #2. (e) Results from all experiments (n = 9, mean ± s.d.) of the time lengths between 
PTZ injections and occurrence of EEG seizure onset, significant gamma power changes and significant hemodynamic changes. 



are needed. A growing body of research indicates that controlling 
seizures may be possible by employing a seizure prevention closed- 
loop treatment strategy'^"". For patients this would be a significant 
breakthrough as they would not be dependent on daily anticon- 
vulsant treatment. Seizure prevention techniques could conceivably 
be coupled with treatment strategies aimed at interrupting the 
process before seizures evolve. Treatment would conceivably occur 
only when needed: on-demand and in advance of an impending 
seizure. 



To develop seizure prevention/intervention models for a specific 
application, it is necessary to first establish that a pre-seizure state 
actually exists. The main finding of our experiments is the early 
changes during the transition from resting interictal state to ictal 
state involves spatiotemporal changes in hemodynamics, as well as 
lEEG gamma power, several minutes in advance of overt seizure 
spikes in local field potentials. Although the mechanisms of seizure 
generation are still unknown, there is evidence that measurable 
changes occur prior to seizures^^'". Accumulated experimental evi- 
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Figure 3 | Comparison between multisite recordings of lEEG signals, gamma power, hemoglobin concentrations and tissue scattering properties, (a) 

Schematic showing the locations of six holes inserted with eight EEC electrodes, (b) Three coronal MRI templates showing eight region of interest located 
on the left and right primary motor (Ml) cortex, hind limb primary somatosensory (SIHL) cortex, somatosensory cortex trunk region (SlTr) and 
anterior thalamus. For each region of interest, lEEG recordings and corresponding gamma powers, hemodynamic responses of [HbT], [Hb02] and 
[HbR] and scattering coefficients (us) responding to the cellular volumes are shown. Results are from rat #2. 



dences have suggested that HFOs in the gamma and fast gamma 
range (40 — 120 Hz) in lEEG recordings, are not only an electrical 
signature of focal epilepsy during seizure- free states, but may also play 
a causal role in the initiation of seizures'. Many neocortical seizures 
were preceded by a build-up of gamma power prior to seizures. On 
the other hand, works on animals and humans indicate that gamma 
band neuronal synchronization in electrophysiological recordings is 
closely related to the BOLD signaP''"*. In parallel to these evidences, 
our findings suggest that DOT can serve as an alternative technique to 
lEEG to detect preictal activities noninvasively. 

Indeed anecdotal reports of preictal changes in brain hemody- 
namics have been reported with SPECT"", fMRI'^, and optical 
imaging^". The dissociating hemodynamic changes from neural activ- 
ities have been discussed in a review*'^ with several hypotheses. One 
hypothesis is that such decoupling is due to the inability of scalp EEG 
or Electrocorticography (ECoG) to detect activities deep in the brain 
or high frequency local field potential activity. According to our 
experiment results, with multisite electrode recordings distributed 



in several regions within the brain, such high frequency changes 
can be observed in some lEEG channels, and came together with 
the hemodynamic changes (Fig. 3b). Our findings of the concurrent 
optical scattering change also support another hypothesis that hemo- 
dynamic responses could be triggered by increased activities of glial 
cells. It has been shown that glial cells play a key role in neurometa- 
bolism and neurovascular coupling''" and cell swelling and shrinkages 
alter the tissue density, thereby changing the scattering properties of 
brain tissues. Thus, glial cell volume changes, which could change 
tissue scattering properties but not detectable in neuronal activity, 
could induce the hemodynamic activities preceding neuronal activity. 
Such decoupling suggests that DOT can be more sensitive than lEEG 
recordings to detect subtle changes. 

Generalized epilepsy is characterized by seizures with diffuse, 
bilateral cerebral involvement, and affects the entire forebrain 
including both hemispheres of the brain, in contrast to the focal 
seizures, which originate from focal brain regions and affect only a 
part of the brain at onsef". In generalized epilepsy multifocal 



SCIENTIFIC REPORTS | 4:3798 | DDI: 1 0. 1 038/srep03798 



5 



resting state 

1 



preictal 

I 



ictal 

1 



ii i i Mi i f ii f li il iiiil i i i i< i i|li l ii iii 



ImI c] n t i i n io 



resting state 

I 



prelctai 

1 




i500nv| mmi ' m ft 



30s 



PP PL 



[!3 lEl EEH 

^ ^ O ET 



3 




resting state preictal 

I 1 



f iU KW iOlllllli 




ictal 

_JL_ 



CI El El 



10 ,10 



0 ■ 0 




Figure 4 | Evolution of [HbT] prior to and during a PTZ induced seizure, (a) shows the positions of five coronal slices, (b), (c) and (d) are the results 
from three seizure experiments (rat #2, #3, and #4) with different seizure intensities. lEEG recordings and statistical activation maps depict [HbT]. 



neocortical regions have the potential to initiate a seizure. Networks 
or neuronal circuits involved in seizures include the seizure focus, the 
initiating circuit, the pathways of spread, and the modulatory cen- 
ters'*^. It has been shown that the interaction between the cortex and 
several nuclei of the thalamus involved in the key circuit for dis- 
charges in the generalized seizures'*^. As can be seen in Fig. 4, during 
early seizure onset, both the thalamus and bilateral cortex are acti- 
vated, and as the seizure progresses, adjacent regions are progres- 
sively activated. 

Despite different patterns of seizures and different origin loca- 
tions, hemodynamic changes all developed from local regions invol- 
ving bilateral cortex and thalamus to the larger area. These findings 
suggest that so-called "generalized" tonic-clonic seizures may be in 



fact localized, at least during the early phase of seizure onset. Indeed, 
the concept of a cortical focus has also been previously discussed in 
other generalized seizure studies'*'"'''^. The concept of focal cortical 
origins during generalized tonic-clonic seizures provides implica- 
tions for understanding and treating this disorder in human patients, 
which may ultimately enable targeted therapies including neurosti- 
mulation, disconnection procedures, and novel pharmacological 
treatments directed at involved regions'*. 

In this study, we did not apply statistical group analysis as has been 
used in other studies^ *. As shown in Fig. 4b-d, for the PTZ-induced 
generalized tonic-clonic seizures model, seizures can occur in a vari- 
ety of lEEG manifestation, even with the same PTZ injection due to 
individual variation for each rat. Such changes in electrophysiology 
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accompany with hemodynamic change with distinct patterns in 
terms of both spatial and temporal distribution. Group analysis could 
potentially smooth or even remove subtle changes in certain areas. 

Heterogeneous patterns of hemodynamic changes have been 
observed in our experiments. Studies of functional activation using 
optical techniques have shown the hemodynamic response in newborn 
infants^^ and a rat modeP*^ to be highly variable. SPECT studies of both 
spontaneous secondarily generalized and electroconvulsive therapy- 
induced tonic-clonic seizures reported heterogeneous patterns of 
changes in CBF during seizures^''^. In fMRI studies^^, a mixed pattern 
of focal increased and decreased BOLD activity has been shown in 
humans and animals during generalized seizures, and relationship 
between neuronal activity and BOLD signals may depend on brain 
region and state. One hypothesis is that preictal vasoconstriction sur- 
rounding the seizure focus could actively shunts blood to the ictal focus 
to prepare the focus for the impending dramatic increase in neuronal 
activity and metabolism^**. However, as shown in Fig. 3, locations with 
either increase (left thalamus and SlTR) or decrease (Ml) of [Hb02] 
all came with significant gamma power increase (t-test, p < 0.05) 
during the preictal period and spikes strong seizure during ictal period. 
Moreover, although it is well understood that the increase in CBF 
caused by the cerebral metabolic rate of oxygen (CMR02) and glucose 
comes with increased neuronal activity, it is still not clear if such CBF 
increase can satisfy local metabolic demands^'^. Depending on the 
interplay between hemodynamics and metabolism, increase of 
[Hb02] could arise when the relative increase in CBF is more than 
the actual need of metabolism. On the other hand, decrease of [Hb02] 
suggests an exceedingly high metabolic demand without appropriate 
compensatory blood flow especially around the seizure focus, rather 
than a decrease in neuronal activity^'*''*^ (Fig. 4b-d). 

While our study provides information for seizure dynamics, it 
does have several limitations that are worth noting and should be 
addressed through future investigation. First, the imaging domain 
should cover the whole brain. Although a large portion of the rat 
brain was included in our current study, we cannot exclude possibly 
important changes in more posterior brain regions. For example, 
much previous work supports an important role for the brainstem 
in generalized seizures both in human patients and animals^'^''*^. 
Second, the studies should be done continuously and over weeks in 
freely moving animals. Also, although the PTZ animal model serves 
the role of a good proof-of-concept model to test our primary hypo- 
thesis, future studies should be focused on chronic epilepsy via an 
animal model with intermittent spontaneous seizures. Third, ana- 
tomical information from MRI should be used in the future study. 
Earlier studies have shown that incorporating subject-specific struc- 
tural priors into the DOT reconstruction process could improve the 
location and quantitative accuracy in human head^^ or breast 
imaging"*^. Besides, subject-specific wearable interface can be built 
with the MRI-derived anatomical information, and applied to freely 
moving rats to facilitate the study of chronic epilepsy with spontan- 
eous seizures. 

In summary, through the use of a PTZ-induced generalized seiz- 
ure model, we were able to investigate early seizure initiation and 
onset as well as seizure propagation patterns. Such studies are essen- 
tial to understanding ictogenesis and guiding seizure prevention 
strategies thereof. To this end, DOT offers an opportunity for observ- 
ing global hemodynamic changes related to neural activities and for 
capturing the genesis and distribution of generalized seizures. The 
possibility of seizure prediction, most importantly, has given hope 
for new warning and therapeutic devices for individuals who cannot 
be successfully treated with current therapies. 

Methods 

All experimental procedures were approved by lACUC committee at the University 
of Florida. 



Animal model. Eleven adult male Sprague-Dawley rats (Harlan Labs, Indianapolis) 
weighing 360-450 g were used for the in vivo experiments. Animals were housed in 
pairs in a controlled environment (12 : 12-h light/dark cycle; food and water ad lib). 

Rats were anaesthetized with i.p. Urethane (1 g/kg, Sigma- Aldrich, Inc.). The 
generalized tonic-clonic seizures were induced by i.p. 130 mg/kg of PTZ (Sigma- 
Aldrich, Inc.)^^. PTZ is the most commonly used GABA^ receptor antagonist used to 
induce absence-like seizures. With low doses (100 mg/kg IP) seizures develop within 
20 minutes after application. PTZ induces all four behavioral phenomena: freezing, 
myoclonic twitches, clonic seizures, and tonic-clonic seizures. Tonic-clonic seizures 
are believed to represent generalized seizures. Although PTZ exerts its action mostly 
via the t-butyl-bicylo-phosphorathionate site of the GABAa receptor^", metabolic 
[14C] 2-deoxyglucose studies demonstrate that there is an uptake or metabolic 
activation in the motor and limbic cortex after PTZ in rat pups and in adult rats, a 
redistribution of glucose uptake from the cortex and hippocampus to the midbrain, 
brainstem, hypothalamus, and septum in adult rats^'. In young rats (<P21), after sub- 
convulsive doses of PTZ, c-fos was expressed in the medial thalamus, cortex, and 
globus pallidus. Occurrence of clonic seizures induced c-fos immunoreactivity in the 
cortex, thalamus, hypothalamus, and brainstem. After clonic seizures, c-fos was found 
also in the striatum, nucleus accumbens, brainstem, and hippocampus*^. 

Time-locked lEEG/DOT were performed continuously for 30 minutes surround- 
ing each PTZ or saline (control experiments) injection, 5 minutes before injection for 
the resting state recording, and then 25 minutes after injection for seizure recording. 
A total of 9/ 1 1 rats were received PTZ, while 2/11 were used as age-matched sham 
controls. 

Electroencephalography. Continuous time-locked lEEG recordings were used to 
confirm the occurrence of seizures. In eight of eleven rats (two rats for control 
experiments; six rats for seizure experiments, including rats #1,#5 — 9 listed in the 
supplementary Fig. Si and Table Si), two electrodes were placed in the bilateral 
frontal cortices (bregma coordinates: AP: 3 mm; ML: ±2.5 mm; DV: 1.5 mm)^^ 
through two holes drilled on the skuU. To monitor and verify the local field potentials, 
multisite lEEGs were recorded in other three rats (rats #2 — 4). In these experiments, 
eight microwire recording electrodes (50 |im polyimide insulated tungsten 
microwire) were implanted through six holes to eight typical bilateral locations on 
Ml, SIHL, SlTr, and anterior thalamic regions: M1(AP: 3 mm; ML: ±2.5 mm; DV: 
1.5 mm), S1HL(AP: 0 mm; ML: ±2.5 mm; DV: 1.5 mm), SlTr(AP: -3 mm; ML: 
±2.5 mm; DV: 1.5 mm), and Thalamus(AP: -3 mm; ML: ±2.5 mm; DV: 5.5 mm). 
Positions of holes and electrodes are shown in Fig. 3. Electrodes for ipsUateral SlTr 
and Thalamus were placed through a same hole since they share same perpendicular 
location. Electrodes were secured using tissue glue. One stainless steel screw electrode 
was implanted as a reference electrode into the midline occipital bone. Local field 
potentials were obtained in differential mode by a neural recording system (RZ5, 
Tucker-Davis Technologies, Inc.) at a sampling rate of 50 kHz and filtered from 0.5 to 
6 kHz. A notch filter was used to remove 60 Hz power line interference. 
Electrocardiogram was also measured to monitor heart rate through one channel of 
this recording system. lEEG data was down sampled from 50 kHz to 600 Hz and then 
band-pass filtered from 2 to 200 Hz. The short-time Fourier transform was 
performed for time- frequency power spectrum density (PSD) analysis using Welch's 
method with a moving window of 6.8 s and 90% overlap. To quantify the level of 
activity in the Gamma frequency (80 — 120 Hz), which is considered closely related 
with the initiation of seizures^, lEEG spectral-band power was calculated by 
integration of power spectrum within that range. To determine the start of significant 
change in gamma power, t-test was performed on each lEEG channel by comparing 
gamma power in 30-second moving window with 80% overlap to the resting state, 
with the nuU hypothesis of none lEEG power change at a significant threshold of p < 
0.05. 

Defining generalized seizure onsets. Typical methods to identify the seizure onsets 
with lEEG are monitoring the appearance of epileptiform lEEG with characteristic of 
high- amplitude or fast activities The start of seizure onset can be determined 
when several number of lEEG channels satisfy such feature. The epileptic seizures 
induced with PTZ in our experiments are generalized tonic-clonic seizures, featured 
with continuous high -amplitude, high-frequency spike and wave activity affecting 
both cerebral hemispheres (covering all lEEG electrodes positions in our experiment) 
from the beginning of the seizure*^. Therefore, to eliminate any confusion regarding 
the onset of the seizure, only the beginning of continuous high- amplitude and high- 
frequency spikes observed in all lEEG channels was considered as the onset of 
generalized seizures, although it could miss some scattered abnormal 
electrophysiological changes in certain channels that occur earlier. Based on this 
standard, seizure onset was identified by visual analysis of electrographic morphology 
and its corresponding PSD for each rat PTZ experiment. 

DOT experiment and image reconstruction. We used a fast multispectral DOT 
system for optical signal recording, as described in detail previously^^. The system is 
schematically illustrated in Fig. 5. Briefly, the system performs continuous -wave 
measurements using 144 light emitting diodes (LEDs) as the excitation sources 
through a time multiplexing approach and highly sensitive photodetectors working in 
parallel for light detection. Forty-eight optical source/detector pairs were distributed 
equally along a cylindrical interface, sequentially shining light to the rat head and 
collecting the diffused scattered light. To ensure good contact and optical coupling 
between the optical fibers/skin, the gap between the head and the fiber interface was 
filled with Intralipid/India ink solution. This DOT system allows a full set of 3-D data 
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Figure 5 | Real time multispectral DOT system for imaging dynamics of generalized seizure. 



collection at 3 NIR wavelengths (660, 780, and 850 nm) for whole brain imaging in 
14.4 Hz. The data acquisition for both DOT and lEEG recording was synchronized 
through a controller board (PCI-7811R, NX). 

In order to provide reliable data for the image reconstruction and analysis, noise 
should be firstly considered and reduced. Basically, instrumental noise, experimental 
errors and systemic physiological artifacts compose the three major sources of noise^^. 
A recent system upgrade to DAQ boards (NI-6358) with higher sampling rate allowed 
an online signal averaging of 100 times for one set of sampling data, which signifi- 
cantly suppressed the random noise as a large part of instrumental noise. A temporal 
low-pass fdter (1 Hz) was applied to attenuate physiological artifacts from move- 
ment, heart beating, and breathing^^-^^. 

After the multispectral signals are collected, an effective model-based reconstruc- 
tion algorithm is critical for accurate image reconstruction in DOT brain imaging. 
Our DOT algorithm is based on the finite element solution to the following photon 
diffusion equation coupled with Robin boundary conditions: 

V-D(r)VO(r) -^i,(r)0(r) ^ -S(r) (1) 

-D(r)Va)n-a<D (2) 

where 0(r) is the photon density, a is a coefficient related to the internal reflection at 
the boundary, /^^(r) are the diffusion and absorption coefficients, and S(r) is the 
source term. The diffusion coefficient can be written as D{r) — + i-i'siy)), 

where //^(r) is the reduced scattering coefficient. The objective of the DOT recon- 
struction algorithm is to recover /i„(r) and /^^(r) at all positions inside the computa- 
tional domain, achieved through a regularized Newton's method to update an initial 
optical property distribution iteratively in order to minimize a weighted sum of the 
squared difference between computed and measured optical data along its boundary. 
A calibration method was applied to reduce the errors caused by the use of different 
source/detection intensities/positions and the system hardware^^. 

Although the continuous wave DOT is theoretically challenged by Arridge and 
Lionheart^^, claiming that simultaneous unique recovery of i.ia{r) and ^.I'Jj) cannot be 
achieved from CW data, our approach is different due to several factors^^. First, 
regularization techniques were used to solve the Hi-posed inverse problem during 
reconstruction. The core procedure in our reconstruction algorithms is to iteratively 
solve the following regularized matrix equation: 

if] + Al)Aq - /^(O*'"^ - O*"^') (3) 

where O'""' and are the measured and computed photon density; /is the Jacobian 
matrix consisting of the derivatives of <!) with respect to D and at the boundary 
measurement sites; I is the identity matrix, and a is the regularization parameter 
determined by combined Marquardt and Tildionov regularization schemes {?. = 0.5 
was used in this study, we have found that the quality of inversion was not sensitive to 



the choice of this parameter). Aq = (ADi,AD2,...,AD;v>A/Vi'^/^«,2'---'A/^rt^)'' is the 
update vector for the optical property profiles, where N is the total number of nodes in 
the finite element mesh. Second, an optimization procedure was applied to search for 
the best parameter set of initial fiaif) and ^(^(r) based on a priori information (bio- 
logically feasible ranges of these initial parameters). The optimization procedure 
seeks best parameter space in the vicinity of the exact one and has allowed the 
minimization of the object function to be confmed in a (/ifl(r),//s(r)) parameter sub- 
space that is close to the exact solution. Using extensive laboratory experiments, we 
have demonstrated that the scattering and absorption coefficients can be distin- 
guished and the solution is unique in many practical scenarios™. These facts have also 
been corroborated by other researchers^^ 

A three-dimensional cylindrical FE mesh, consisting of 22,764 tetrahedron ele- 
ments, was used for the image reconstruction. To efficiently process a large amount of 
data containing thousands of time points for one set of experiment, a GPU (graphic 
processing units) based parallel code^^ was implemented using CUDA programming 
model (NVIDIA), allowing the time consuming image reconstructions for one set of 
experiment to be processed within one day. 

DOT analysis and rat atlas registration. Assuming oxy- and deoxy- hemoglobin 
(Hb02 and Hb^) are the two major chromophores in the rat brain, the absorption 
coefficients at 660 nm, 780 nm and 850 nm were used to obtain [Hb02] and [HbR] 
using the Beer-Lambert law^^ and least square fitting through pseudo inverse matrix 
calculation. The concentration of total hemoglobin (HbT) was then calculated by 
summing up [Hb02] and [Hb^] which is proportional to CBV. [HbT] has been 
suggested to map cerebral activity in NIRS for its better spatial specificity than [Hb02] 
or [HbR]^*. Activation maps were generated by tracking the change of [HbT] that 
exceeded significant threshold, because HbT contains the information of both 
[Hb02] and [Hh^] and reflects the local CBV. For statistical analysis, thirty seconds 
(432 images) consecutive [HbT] images preceding PTZ injection were chosen as a 
baseline period to which all analysis periods were compared. Several time series 
images, each covering a 10 seconds time bin (144 images), were selected from the 
interictal, pre-seizure, and seizure periods. A t-test was performed by comparing each 
time bin to baseline using a voxel wise 2-sample t-test, with a significant threshold of p 
< 0.05. Activation maps were then overlaid with template MRI images obtained from 
the Karolinska rat atlas (http://expmr.ki.se/research/ratatlas.jsp) for anatomical 
guidance. 

Activation maps were visualized through the following procedures. First, Amira 
(Visage Imaging, Inc.) was used to detect the brain surface from the MRI template and 
generated a brain mask. Then, DOT image coordinate was aligned with the atlas 
coordinate by comparing the measured bregma position from experiments and its 
according position in the atlas. Next, DOT images within the brain mask were 
interpolated to each voxels (1 mm^) in the atlas using Matlab. Results were converted 
and stored in NIFTI format, which were visualized with MRJcron (http://www. 
mricro.com/). 
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